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The identification of HCN and HNC in Carbon Stars: 
Model Atmospheres, Synthetic Spectra and Fits to 
Observations in the 2.7-4.0 /im Region 
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ABSTRACT 

Model carbon star atmospheres and synthetic spectra have been c alculated using 
the recent HCN/HNC vibration rotation linelist of lHarris et"ai"ll2002bl ApJ, 578, 657. 
The calculations are repeated using only HCN lines and show that HNC has a signif- 
icant effect upon the temperature, density and optical depth of a stellar atmosphere. 

We fit synthetic spectra in the 2.7-4.0 /itm re gion to observed ISO spectra of 
the carbon stars WZ Cas and TX Psc obtained bv lAoki et aTlll998t A&A, 340, 222. 
These fits allow us to identify absorption by HNC in the spectrum of WZ Cas at 
2.8-2.9 /mi, and to determine new independent estimates of effective temperature and 
log(iVc/iVo) ■ The findings reported here indicate that absorption by both HCN and 
HNC is needed to fully explain the observed stellar spectra and represent the first 
identification of HNC in a star. Q branch absorption by the HCN A«2 = 1, A«3 = 1 
and Ai>i = 1, Av2 = — 1 bands at 3.55 and 3.86 /im respectively, are identified in the 
spectrum of WZ Cas. 

Key words: StarsAGB, stars: carbon, stars: atmospheres, infrared: stars, molecular 
data. 



1 INTRODUCTION. 

It is well known that HCN is an important opacity source 
in carbon stars and one which can dramati cally affect the 
struc t ure of the atmosphere s of these stars jEriksson et alJ 
Il984l : Ijflrgensen et alJll985l) . It is therefore vital that the 
opacity of this particular molecule should be accurately ac- 
counted for. Several recent works have used different sources 
of HCN opacity data. The recent carbon rich model atmo- 
sphere calculations of lAoki. Tsuii fc Ohnakal (Il998l I1999T) 
used the HCN band transition dipol es and molecular con- 
stants dete rmined expe r iment ally by iMaki. Quapp fc Kled 
lll995l) and IMaki et alJ (Il996l) . to accou nt for HCN opac- 
ity. In contrast other recent works llLoidl et alJ Il999t 
[ Loidl. Lancon fc Jorgensenl 1200 it iJor^ense^i^ton^^jOidjl 
l20^oF have used the ab initio data of ljoTgeirseri^t^lTiH^S^^ 
Although ab initio data cannot rival laboratory data in ac- 
curacy, laboratory data are limited to tiny fraction of the 
billions of possible HCN transitions. As a consequence the 
laboratory data cannot possibly account for all of the HCN 
opacity. This makes ab initio data more useful for high tem- 
perature calculations which require a more complete descrip- 
tion of molecular vibration rotation spectra and hence opac- 



ity. The ab initio data of |j0rgensen et al.l (Il985l) was the 
first extensive HCN opacity data set, and was a step for- 
ward from the empir ical estimates of HCN opacity used by 
lEriksson et alJ (1984). However over the intervening decade 
and a half the advances in molecular quantum mechanical 
techniques coupled with increasingly more powerful comput- 
ers has improved the quality and quantity of spectroscopic 
calculations on triatomic molecules. These improvements en- 
a bled the calculation of the extens i ve HCN and HNC linelist 
of lHarris. Polvanskv fc Tennyson] (l2002tJ) . which we use in 
this work. The linelist of Trhu^^t£dTT2002b[) is both more 
extens ive and accurate than the data of lj0rgensen et alJ 
(1985). 



A further adva ntage of the linelist of |l arri^^^l 



( 2002b]) over those of lJargensen et all £l985) and Aoki et al 
(1998) is that it also takes into account the effect 
of the isomer of HCN: HNC. Recent calculations by 
iBarber. Harris fc Tennvsonl §002) have yielded a HNC to 
HCN ratio, at thermodynamic equilibrium, of 0.13 at 2600 
K and 0.19 at 3000 K. This coupled with the fact that the 
band transition dipoles of the fundamental bands of HNC 
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are stronger than their HCN counterparts, imply that HNC 
will contribute to the opacity of C-star atmospheres. 

To determine the effect the new opacity data of 
lHarris et alJ (l2002bl) has on model C-star atmospheres and 
synthetic spectra, we present a grid of model atmospheres 
calculated using this data. We then discuss these model at- 
mospheres and compare them with existing C-star model 
atmospheres. We compare the synthetic spectra with exist- 
ing ISO observations of the C-stars TX Psc and WZ Cas in 
the 2.7 to 4 /im region. This spectral region lies between the 
two strong CO bands, Av=l and 2, and contains several 
discernible HCN bands and the HNC H-N stretch funda- 
mental band. From this comparison we determine effective 
temperatures and C/O ratios for these star. The effect of 
HNC opacity is identified for the first time in the spectrum 
of WZ Cas, and we attempt to identify Q branches of two 
HCN modes in the spectra of WZ Cas. 



2 HCN AND HNC. 

The HCN/HNC linelist of lHarris et all (|2002bh covers tran- 
sitions between all HCN and HNC energy levels below 18 000 
cm -1 above the HCN zero point energy and with angular 
momentum quantum number J sC 60. There are approach- 
ing 400 million lines and around 125 000 vibration rotation 
energy levels in this linelist. The accuracy of the results of 
this calcu lation room tempera ture situations has been dis- 
cussed in lHarris et all (l2002al) . In particular it was found 
that the band transition dipoles calculated from this linelist 
were within combined systematic and experimental error. 
All the fundamental band centres are displaced from the 
laboratory values by less than 4 cm -1 , except the HNC H-N 
stretch fundamental which has an error of 12.5 cm -1 . The 
rotational error, corrected for vibrational error, for the lines 
of the fundamental vibra tional modes of HCN/HN C are less 
than 0.7 cm" 1 at J=20 llvan Mourik et al-lEoOlL Many of 
the lines in the linelist are extremely weak and make no 
real contribution to the overall opacity. We therefore reduce 
the number of lines used in the calculation of the model at- 
mospheres and synthetic spectra by cutting out lines based 
on absolute intensity. We used a minimum absolute inten- 
sity at 3000 K of 10 -15 cm _2 atm _1 , the strongest lines are 
of the order of 0.1 cm _2 atm _1 . This truncation dramati- 
cally reduced the number of lines, with little effect on the 
overall opacity. For the model atmospheres lines were also 
selected on the basis of proximity to one of the frequency 
points used in the calculation. We must correc t an error in 
the valu e of the constant given in equation 1 of lHarris et alJ 
(2002b). This constant should read 3.56589 xlO 7 , and not 
3.50729xl0 -6 . 

To clearly show the effect of HNC opacity on the atmo- 
spheres of C-stars it i s necessary to separa te the combined 
HCN/HNC linelist of lHarris et all (l2002bl) into HCN and 
HNC data. There ar e 125 000 energy levels in the linelist of 
lHarris et al ] j2002bTl . it would be very time consuming and 
tedious to identify by hand energy levels as HCN or HNC. 
Consequently we used the identificatio n algorithm used by 
lHarris. Polvanskv fc Tennvsonl i2002aT) . In this method in- 
dividual energy levels are labelled HCN or HNC on the ba- 
sis of the strengths of the transition from that level to the 
vibrational ground HCN and HNC states. An energy level 



that has a transition dipole to the ground state of HCN 
which is several orders of magnitude stronger than the tran- 
sition to the ground HNC state, can confidently be identi- 
fied as HCN and vice-versa. Two problems arise with this 
method. The first and most important involves energy lev- 
els which have very weak transitions to the ground state. 
For example the HCN (0,8,8,0)^(0,0,0,0) band is forbidden 
by the normal dipole transition rules for a linear triatomic 
(Al = —1, 0, +1). This transition is so weak that it is below 
the noise level in the calculation, states with a transition 
to the ground state this weak we do not assign as HCN or 
HNC. Secondly, the idea that HCN and HNC are separate 
molecules is incorrect at high energies. The wave functions 
of some high lying states have magnitude in both HCN and 
HNC potential wells and so cannot be described as HCN or 
HNC, these are called 'delocalised' states. Such states can 
have transitions to the HCN and HNC ground states which 
are of roughly the same order of magnitude. 

The calculations that we make in this work which are 
described as 'HCN only', exclude only the HNC levels which 
we can confidently identify. There may be other HNC energy 
levels in the data which are unassigned. However this 'HCN 
only' data will exclude the strongest HNC transitions and 
will allow us to identify the effect of HNC on our model at- 
mospheres and spectra. We label energy levels confidently 
identified as HCN as i~0, HNC levels as 1=1, delocalised 
states as i = 2 and unassigned states as i — 3. Figure 
shows a plot of HCN and HNC absorbance at T=3000 K 
as a function of wavelength, from 2.6 to 4.1 /im. The three 
curves show absorbance calculated with all lines, absorbance 
calculated with HCN lines only and absorbance calculated 
only with confidently identified HNC lines. Over this spec- 
tral range HNC absorption only exceeds that of HCN at 
wavelengths between 2.7 to 2.9 /im. However between 3.2 
and 3.9 /im HNC appears to contribute significantly (10- 
15%) to the total opacity. The effect of this HNC absorption 
can be seen in our synthetic spectra, see sections 15.2.11 and 

EH 

We ha ve used the recent ab initio HCN/HNC partition 
function of Ba rber et al.l (120021) in all of our calculations. 
This partition function was calculated by direct summation 
using the HCN/HNC energy levels of lHarris et"aT] j20O2b) 
up to J=32 and includes an extrapolation of these levels 
up to J=91, which is enough to rotationally converge the 
partition function. The error of this partition function has 
been estim ated to be 0.6% at T=2600 K and of 0.8% at 
T=3000 K jBarber et alJl2002l) . HNC energy levels are also 
included in the calculation. 



3 STELLAR SPECTRA. 

The spectral region required for this program is difficult to 
obtain from ground-based observatories and is well matched 
to ISO capabilities. In particular the SWS instrument on 
ISO provides sufficient resolution to resolve individual vi- 
bration rotation bands. The observations used for this study 
were taken as part of the Japane se guaranteed ob serving 
time (REDSTAR1, PI T.Tsuji, e.g.. lAoki et~ai1 dl998h l. This 
programme obtained a number of infrared measurements of 
Carbon stars. For the purposes of our study we selected the 
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Figure 1. The wavelength dependent absorbance of HCN and HNC. The three curves represent the absorbance of HCN/HNC (all lines), 
the absorbance of HCN (all lines except (i' = 1, i" = 1) and HNC (lines with %' = 1, i" = 1). The main HCN and HNC features are 
labelled. See text for a full description of i. 



high resolution observations made for TX Psc (N-type) and 
WZ Cas (SC-type) from the ISO archive facility. 

The observations were made in with the Short Wave- 
length Spectrometer (SWS) in its full resolution grating 
mode known as SWS06. Here we use spectra with a res- 
olution of 1000-2000 from bands IB, ID and IE covering 
the wavelength ranges 2.74-3.02, 3.02-3.52 and 3.52-4.08 
pirn. The observational data are scientifically validated and 
were processed using the standard archive reduction soft- 
ware. The spectra were inspected using the ISO Spectral 
Analysis P a ckage to check for glitches and bad data points. 
lAoki et alJ lll99fet) note that the flux discrepancies are at 
most 10% for the REDSTAR1 dataset. For the subset of 
data used for this paper we find the flux calibration errors 
between SWS bands are at the level of a few percent, as 
expected for band 1 SWS observations. The wavelength cal- 
ibration for the SWS instr ument is measured t o be within 
1/8 of a resolution element ISalama et alJ ll2000|) . 



4 PROCEDURE 

4.1 Ionisation-dissociation equilibrium 

In atmospheres of cool (T e ff < 6500 K) carbon rich stars 
the ionisation-dissociation equilibrium is governed by the 
CO molecule, that is, the densities of carbon-containing 
species depend on the C/O ratio (see Tsuji 1973 for more 
details). The ionisation-dissociation equilibrium is computed 
for ~ 100 of species including the most abundant carbon- 
containing molecules. Const ants of chem ical equilibrium for 
most species are taken from iTsuiH lll973t) and constants for 
HCN are converted into o ur model atmospher e format using 
partition function data of lBarber et alJ l2002f) following the 



scheme of lPavlenkol <2002bl) . HCN and HNC molecules are 
treated as a single molecule with a single partition function 
and dissociation potential. 



4.2 Opacities 

The stand ard set of con tinuum opacities included in AT- 
LAS12 bv lKuruczl lll999l) are used. Computations of a few 
opacities were add ed or u pgraded: the absorption of C~ 
llMverscoug hll968t) and H^ <Dovlell968ft . To compute opac- 
ities due to b ound-free absorption of C, N, O atoms we used 
TOPBASE iSeaton. Zeipen fc Tullvl Il992l) . (cross-sections 
are available on http://www .mao.kiev.ua/staff/yp/), see 
JPavlenko fc ZhukovskJ 2003) for more details. 

To compute opacity due to bound-bound absorption of 
atoms and molecules in atmospheres of C-giants we used 
linelists from different sources: 

(i) Electronic and vibration rotation line data of diatomic 
molecules (CN C 2 , MgH, H 2 , SiH, CO) were taken from 
CDROM 18 of iKuruczj <ll993l) . For synthetic spectra com- 
putations we use d the vibration rotation lines of CO from 
iGoorvitcbl \l994) . 

(ii) Atomic lin e data were taken from VALD 
ftCupkaet alJll999l) . 

(iii) HCN / HNC v ibration rotation lines were taken from 
lHarris et all d2002bl) . see section H 

(iv) CS vibration rotation lines were taken from 
IChandra et all lll995l) . 

Other works <Aoki et alJ Il998t Il999t IJorgensen et alJ 

2000), have included C2H2 opacity. For warmer C-stars, 
such as WZ C as and TX Psc, C2H2 has only a small ef- 
fect on opacity llAoki et al.lll998l) . Hence it is not necessary 
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to include H2C2 within our calculations. However for cooler 
C-stars with high log(iVc/-/Vo) , H2C2 may well be more 
important and should be included within the calculation. 

For the model atmosph ere calculations we used 
the op acity sampling approach ISneden. Johnson fc Kruprl 
( 1976) to account for the absorption by atomic and molec- 
ular lines. Absorption of molecular electronic bands: CaO 
(C X E - X^), CS(A 1 E - X 1 ^), SO (A 3 n - X 3 E), SiO (E X E- 

x j e), siO(A 1 n - x 1 s+, no (c 2 n r - x 2 n r ), NO(B 2 n r - 

X 2 n r ), NO(A 2 E+ - X 2 n r ), MgO(B 1 E+ - X 1 E+), Aio(c 2 n 
-X 2 E), A10(B 2 E+ - X 2 E+) are taken into account by the 
just overlappin g line approach (JOLA), using the subpro- 
grams of Nersisian, Sh avrina fc Ya kcmchuk (1989]). In gen- 
eral, these molecules absorb in the blue part of the spectrum 
and their bands are not principal opacity sources in stellar 
atmospheres. Nevertheless, the JOLA algorithm enables us 
to reproduce both the relative strength and sha pes of the 
spectral en ergy distributions for late type stars jPavlenkd 
ll997l.l2000D . 

4.3 Model atmospheres 

We computed a grid o f model atmospheres of C -giants with 
the SAM12 program iPavlenkd J2002al. 12003ft which is a 
modification of ATLAS 12 (Kurucz 1999). SAM12 is a clas- 
sical model atmospheres code: plane-parallel media, LTE, 
no sources and sinks of energy within the atmosphere, en- 
ergy transferred by radiative field and convection. Convec- 
tion is treated in the frame of mixing length theory, in our 
computations we adopt l/H = 2.0, convective overshooting 
jKurucdll999l) was not co nsidered. Opacity sampling treat- 
ment Jsned^nat alJll976T) is used to account for atomic and 
molecular absorption. We investigate the sensitivity of our 
synthetic spectra to changes in effective temperatures T c fr , 
abundances of carbon log N(C) and mi croturbulent veloci- 
ties V t. We adopt the "solar values" of I Anders fc Grevessd 
( 1989) for elements other than carbon, for example oxygen 
log N(O) = -3.12 is adopted where J2 N > = L 

Model atmospheres were computed for effective tem- 
peratures (Tog = 2500 - 3500 K), carbon abun- 
dances log N(C) = -3.19 to -3.09, carbon isotopic 
ratio 12 C/ 13 C = 3 -10. Around 100 model atmo- 
spheres were generated for this work. The models com- 
puted with the complete HCN/HNC linelist are avail- 
able from http://www.mao.kiev.ua/staff/yp/Results as 
files Mod.CglO.tar.gz and Mod.Cg03.tar.gz. 



4.4 Synthetic Spectra. 

Numerical computations of theoretical radiative fluxes F v 
were carried out within the same classical appro ach as our 
mode l atmospheres using the WITA6 program jPavlenkd 
2000). The shape of every line was determined using the 
Voigt function H(a,v). We account for natural broadening 
as well as van der Waals and Stark broadening of lines. Due 
to the low temperature and densities of C-giant model atmo- 
spheres Stark broadening may be ignored. Radiative broad- 
ening may play a major role in the outermost layers; in deep 
photospheric layers van der Waals broadening becomes im- 
portant. Dampi ng constants were t aken from line databases 
such as VALD jKupkaet al.l [t999) or computed following 



lUnsoldl Jl955T> . The generation of our synthetic spectra are 
performed using the same opacity source lists and sets of 
abundances, as our model atmospheres. This approach al- 
lows us to directly see the impact of abundance changes on 
the temperature structure of the model atmosphere as well 
as the synthetic spectra. Computations of synthetic spectra 
were carried out with wavelength steps of 0.05 nm. 

4.5 Fit to observed spectra 

There are several parameters that can be adjusted to fit syn- 
thetic spectra to observed stellar spectra. The most impor- 
tant of these include the logarithm of the carbon to oxygen 
ratio (\og(Nc/No) ), effective temperature (T of f), surface 
gravity (logg), metallicity (\og(Z / Zq)) and turbulent veloc- 
ity (Vt). Fits obtained by varying these paramete rs are of- 
ten not unique, for example lJorgensen et al.l i2000f) obtained 
two good fits for TX Psc, with different values of these pa- 
rameters. In this work we attempt to gain good fits to the 
observed stellar spectra with \og(Nc/No) and T e g , but if 
necessary we change Vt and the metallicity. 

To obt ain a fit we fix lo g.q to 0.0, this value was es- 
timated by lAoki et alJ lll998l) from a measured radius of 
200 Rq for TX Psc and an assumed mass of between 1 
and 2 Mq, We also fix Vt to 3.5 km s _1 and metallic- 
ity to solar values, then calculate synthetic spectra on a 
log(iVc / 'No) , Tog grid. When a reasonable fit is found we 
adjust log(A r c/A r o) and T G g to obtain a better fit and if 
necessary change Vt and metallicity to further improve the 
fit. We did not find it necessary to change log g to obtain 
good fits. Normalisation of the synthetic spectra and stellar 
spectra was carried out by dividing through by the Planck 
function of the relevant T c ff , then scaling the spectra so that 
the maximum flux is 1. The temperature of the black body 
spectrum used to normalise the stellar spectra is the T e fi of 
the best fit synthetic spectrum. The division of the observed 
spectra by the Planck function is adequate to normalise the 
spectrum to an pseudo-continuum. 



5 RESULTS 

The calculations re ported here are the first to use the new 
HCN/ HNC data of Harris et alJ J2002bl) . The fact that the 
data of lHarris et alJ ( l2002bl) contains more lines implies a 
greater atmospher ic opacity than accoun t ed for in precedin g 
HCN/HNC data Jjargensen et alJll98Et lAoki et al]ll99Sft . 
This increased opacity may have a significant effect on both 
the structure and spectra of our model atmospheres and 
synthetic spectra. 

5.1 Model Atmospheres 

Figure |2] shows a plot of pressure versus temperature for 
9 different model atmospheres with a log(iVc/iVo) of 
0.01, temperatures of 2600, 2800 and 3000 K, and with 
HCN/HNC opacity with only HCN opacity and without 
HCN/HNC opacity. Symbols are placed on the curves at 
value of log tr of -4, -3, -2, -1, 0, 1 and 2. The difference be- 
tween the model atmospheres with and without HCN/HNC 
opacity at 2800 and 2600 K is considerable. At a given op- 
tical depth the atmospheres with HCN/HNC opacity are 
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cooler and at a much lower pressu re, which is consisten t 
with the trends originally reported bv lEriksson et al.Hl984) . 
However in the model atmospheres with an effective tem- 
perature of 3000 K, most of the HCN/HNC has been de- 
stroyed and the differences between spectra with and with- 
out HCN/HNC opacity is much reduced. None-the-less these 
differences are still significant, especially in the outer parts 
of the atmosphere. The model atmospheres with HCN/HNC 
and with only HCN, have discernible differences. At a given 
optical depth the models with HCN opacity only are hotter 
and at higher pressures than the models with HCN/HNC 
opacity. Much of this change in the optical depth pressure- 
temperature relationship is due to the increase in opacity 
and hence optical depth as a result of HNC absorption. 
However some of the differences are also attributable to 
changes in the pressure temperature curve of the atmo- 
spheres, brought about by increased absorption. The change 
in the pressure temperature curve is most pronounced in the 
2600 K models. 

Figure [3] shows a plot of pressure versus temperature 
for 9 model atmospheres with \og(N c /N ) of 0.003, 0.01 
and 0.03 at a temperature of 2800 K, with and without 
HCN/HNC opacity and with HCN opacity only. As in figure 
|2]the model atmospheres with HCN/HNC opacity are cooler 
and at much lower pressure than those without HCN/HNC 
opacity. It is also clear that the atmospheres with a higher 
C/O ratio have a hotter low er pressure atmosp here, this is 
consistent with the findings of lAoki et alJ (ll998T> . The reason 
is that more carbon rich molecules can form at higher C/O 
ratio resulting in a more opaque, so hotter and less dense, 
atmosphere. All these pressure-temperature cur ves are con- 
sistent wit h the model atmosphere s reported bv lAoki et alJ 
l|l99Sl) and IJorgensen et alJ feOOOl) . A gain there are differ- 
ences between the models which account for HCN/HNC 
opacity and the models which account for only HCN opacity. 

5.2 Spectra 

Figure^Jshows synthetic spectra calculated at effective tem- 
peratures of 2600, 2800, 3000, 3100 K. At T cff of 2600 and 
2800 K the atmospheric temperatures and densities are con- 
ducive to the formation of HCN/HNC, so the spectra be- 
tween 2.7 and 4.0 /jm are dominated by HCN absorption. 
Between 2.9 to 3.4 /jm is the characteristic strong absorp- 
tion feature which is predominately caused by the HCN, 
H-C stretching mode. There are strong narrow absorption 
features around 3.6 /im caused by HCN Q branches, see 
below. Immediately either side of the 3.6 /im feature are 
broader weaker absorption features corresponding to P and 
R branches of the same bands. There are also strong, nar- 
row absorption features longwards of 3.8 /jm, caused by a 
combination of CS and HNC bands, see below. Finally the 
broad absorption around 2.8 /im is caused by the HNC, H-N 
stretching fundamental and its hot bands. At T c g = 3000 K, 
all these absorption features are weaker, the narrow features 
at 3.6 /im and longwards of 3.8 /im have all but vanished. 
The 3 to 3.4 /im feature has narrowed and weakened. The 
HNC absorption at 2.8 /im is considerably weaker, and the 
spectrum in this region is becoming dominated by CO. At 
T c ff=3200 K, almost no trace of HCN and HNC absorption 
is left, even at 3 /im. Longwards of 2.9 /im the spectrum is 
flat. These changes are a direct result of the temperature de- 



pendence of the HCN and HNC abundance. As the effective 
temperature increases to 3000 K the HCN/HNC molecules 
begin dissociating to form CN, this process is complete by 
3200 K at which point little sign of HCN is present. 

Figure shows synthetic spectra calculated with an ef- 
fective temperature of 3000 K and with log(iVc /No) of 0.03, 
0.02, 0.01, 0.003 and 0.001. At T cff =3000 K the HCN/HNC 
abundance is highly sensitive to both the C/O ratio and 
temperature. This is seen in both figures 01 and |5] The 3-3.4 
/xm feature is strong for \og(Nc/No) between 0.001 and 
0.01, but is weak at log(A r c/A r o) = 0.02 and almost no ex- 
istent at log(Nc/No) =0.03. Here higher C/O ratios give 
rise to a greater abundance of carbon rich molecules which 
in turn makes the atmosphere more opaque and so hotter 
and less dense, see figure |3 Such conditions do not favour 
the formation of HCN/HNC. So a high C/O ratio can result 
in the destruction of HCN/HNC. These sensitivities to the 
C/O ratio at T e ff~3000 K, make it difficult to fit spectra to 
the higher temperature C-stars, such as TX Psc, see section 

EH 

As discussed in section|3]we have obtained suitable spec- 
tra of the carbon rich AGB stars WZ Cas and TX Psc, 
from the ISO archi ves. The inclusion of the new extensive 
HCN/HNC linelist jHarris et all2002bl) within our synthetic 
spectra allows the identification of HCN and HNC features, 
within these stellar spectra. The new HCN/HNC are more 
complete and in the case of ab initio data more accurate, 
than the preceding HCN/HNC data and so allow the deter- 
mination of new effective temperatures and C/O ratio for 
these stars. 

5.2.1 WZ Cas. 

Figure HJ shows the spectrum of WZ Cas together with the 
best fit synthetic spectrum calculated with HCN/HNC opac- 
ity and a synthetic spectrum calculated using only HCN 
opacity. Both of the synthetic spectra were computed for a 
model atmosphere of T e fi =2800 K, a turbulent velocity of 
3.5 km s" 1 and a value of \og(N c /N ) of 0.003. To permit 
easy identification of the broad absorption features within 
the spectra, both synthetic and observed spectra are shown 
at a resolution of 1.5 nm. Recent measurements and fits for 
the effective temperature and log(A r c/A r o) are shown in 
table □ There are large differe nces between the values of 
Teg. The estimates of T e e from lDvck. van Belle fc Benson! 
( 1996) is 350 K hotter than the estimate of T e ff from our 
models. However if we use this T e s in our models, HCN 
does not form in large enough quantities to form strong 
absorption features, see figure [I] There are large e r rors o f 
200 K, on the value of T cff quoted bv lDvck et aD il99fil> . 
Our estimate o f the T e g and log(N c/No) for WZ Cas is 
closer to that o flLambert et all <ll986T) and that of lAoki et alJ 
( 1998) whi ch is an recalibration of the infrar ed flux method 
fl RFM) iBlackwell. Petfqrd fc Shaflis! ll98fJl determination 
bv lOhnaka fc Tsuiil 

In general the major features of the synthetic spectrum 
with both HCN and HNC opacity, fit well the observed spec- 
trum of WZ Cas. The absorption feature above 3.8 microns, 
is attributable to a combination of the Q branch absorption 
of the (0,1,1,0)— >(1, 0,0,0) band together with its hot bands 
and the Av — 2 bands of CS. The general shape of this fea- 
ture is well reproduced in both synthetic spectra. However 
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Figure 2. Plot of pressure versus temperature on logarithmic scales for models with HCN/HNC lines (solid, + symbols), with HCN 
lines only (dotted, X symbols) and without any HCN/HNC lines (dashed, o symbols), at T=2600 (thick), 2800 (medium) and 3000 K 
(thin), with log(Nc/No) of 0.01. Symbols are placed at Rosseland optical depths of logT/j=-4, -3, -2, -1, 0, 1, 2, 3. 




Figure 3. Plot of pressure versus temperature on logarithmic scales for models with HCN/HNC lines (solid lines, + symbols), with 
HCN lines only (dotted lines, X symbols) and without any HCN/HNC lines (dashed, o symbols), at T=2800, with \og(N c /N ) of 0.003 
(thick), 0.01 (medium) and 0.03 (thin). Symbols are placed at Rosseland optical depths of logTR=-4, -3, -2, -1, 0, 1, 2, 3. 



the positions of the individual Q branches in the synthetic 
spectra are slightly displaced from the same Q branches in 
the observed spectrum. This is a resu lt of the errors on th e 
line frequencies in the ab initio data of lHarris et alJ ((2002b). 
There is slightly weaker absorption in this region in the syn- 



thetic spectrum with only HCN absorption than there is in 
the spectrum which uses the full HCN/HNC opacity. This is 
most likely a result of the effect of absorption by the HNC 
(0,0,0,0)— >(0, 1,1,1) band, centred on 4.03 /im, and its hot 
bands. From figure^it is expected that HNC will contribute 
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log(N c /N o )=0.01 
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Figure 4. Synthetic spectra with a log(AT c /jV ) of 0.01 and effective temperatures of 2600, 2800, 3000 and 3100 K. 



T=3000 K 




2.8 3 3.2 3.4 3.6 3.8 
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Figure 5. Synthetic spectra with an effective temperature of 3000 K and log(N c /N ) of 0.03, 0.02, 0.01, 0.003 and 0.001. A constant 
has been added to the flux of some of the spectra, to stagger them. 



Table 1. Effective temperatures, C/O ratios and surface gravities for WZ Cas. 



Work T cff log(7V c /Ar ) log(g) Notes 

Lambert et al. ( 1986) 2850 0.0043 T e g from colour index's, \og(N c /N ) from models. 

Dvck et al. (19961 3140±193 interferometry 

Aoki et al. (1998') 2900 0.0043 0.0 Recalibrated IRFM determination of Ohnaka fc Tsuii ( 19961 

This work 2800 0.003 0.0 model 
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2.8 3 3.2 3.4 3.6 3.8 4 

Wavelength (microns) 

Figure 6. Synthetic spectra of carbon rich atmospheres at T e g =2800 K, logf Ng /No) of 0.00 3, with a Vt=3.5 km s _1 in the region 
2.7-4.0 ^m. The spectra were calculated with the full HCN/HNC linelist of IHarris etahl <2002ll and with only HCN lines selected. For 
comparison the observed spectrum of WZ Cas is shown. The resolution of these spectra is 1.6 nm. 



to the overall HCN/HNC opacity in this region. Some of the 
changes in the synthetic spectra may also have been brought 
about by changes in the optical depth, pressure and temper- 
ature of the model atmosphere. 

The absorption in the 3.6 /im region is attributable to 
the Q branches of the HCN (0,0,0,0)^(0,1,1,1) band. Again 
the general shape of this feature in both synthetic spectra is 
well reproduced, but the Q branches in the synthetic spectra 
are displaced from the observed spect ra, a result of t he erro r 
on the ab initio HCN/HNC data of IHarris et"aD ]2002bl) . 
Shortwards of 3.55 /im the spectrum which contains both 
HCN and HNC absorption agrees well with the observed 
spectra. The spectrum which does not contain HNC data 
appears to be missing absorption, especially in the 2.9 fim 
region. The broad absorption feature from 2.95 to 3.4 /zm 
which is widely attributed to the HCN H-C stretching fun- 
damental and hot bands, centred on 3.02 /mi, and its hot 
bands. This region has bee n well documented in the past (for 
example lAoki et al] lll99fet0 . There should also be a little ab- 
sorption from the C2H2 symmetric and anti-symmetric H-C 
stretch modes. However the difference between the intensity 
of this feature in the synthetic spectra with and without 
HNC, suggest that much absorption in this region is a re- 
sult of HNC. The strongest HNC mode in this region is the 
(0,1,1,0)— ►(1,0,0,0) band, centred on 3.125 /tm, and its hot 
bands. The effect of this band on the combined HCN/HNC 
opacity can be seen in figure 

The biggest difference between the synthetic spectra 
with HCN/HNC lines and with HCN lines only, lies in 
the 2.85 /xm region. The dominant absorption in this re- 
gion of the synthetic spectrum with only HCN opacity, is 
caused by CO absorption from the Av = 2 band. How- 
ever, in the synthetic spectrum which includes the full 



HCN/HNC opacity there is additional absorption. This is 
the HNC (0,0,0,0)^(1,0,0,0) fundamental and hot bands, 
centred on 2.73 /jm. From the ab initio absorption spectrum 
of HCN/HNC (figure 01 it would be expected that if HCN 
and HNC are in thermodynamic equilibrium there should be 
significant HNC absorption in this region. In the 2.85 /im re- 
gion the difference between the observed and synthetic spec- 
tra are also at their greatest. The synthetic spectrum which 
includes the full HCN/HNC linelist fits both the shape and 
flux of the spectrum of WZ Cas considerably better than 
does the synthetic spectrum which contains only HCN ab- 
sorption. Thus it is here that we identify the effect of HNC 
absorption in the spectrum of a carbon star for the first 
time. 

The synthetic spectra of lAoki et al] lll998l) were cal- 
culate d using available laboratory HCN data. lAoki et al] 
l|l998h were not able to fully account for the HCN absorp- 
tion between 3.8 to 4.2 pm, using only this labor atory data. 
To gai n a closer agreement with observed spectra FAoki et ail 
( 1998) had to empirically estimate opacity from Avi = 2, 
Av2 = — 1 hot bands with V2 > 5. Our synthetic spectra do 
not suffer from this data deficit and are on the whole are in 
good agreement in this region. The minim um frequency of 
the synthetic spectra of lAoki et al.l ijlOOST) is 2.9 /im, so we 
are unable to comment upon the spectra at 2.8 /im where 
HNC is a strong absorber. 

For the spectral ranges 3.55 to 3.63 /zm and 3.83 to 
3.97 /mi we plot the observed spectrum of WZ Cas together 
with the best fit synthetic spectrum in figures and [£] re- 
spectively. To help the reader identify the bands of HCN 
responsible for the absorption features in the synthetic and 
hence observed spectra we have plotted the ab initio HCN Q 
branch lines at the bottom of these figures. The band centres 
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of the ab initio data of lHarris et alJ J2002bl) are identified 
by the lower marks and digits and the experimentally de- 
termined band centres by the upper marks and digits. The 
indexes of the band centres are given in tables|5|and tabled 
To show the relatively broad absorption features the resolu- 
tion of both the observed and synthetic spectra has been set 
to 0.75 nm. The typical spacing between Q branch lines in 
triatomics is considerably smaller (about 0.1 nm) than the 
spacing between P and R branch lines (about 4 nm). Conse- 
quently, the intensity of the P and R branches is distributed 
over a large wavelength range. Hence the P and R branches 
are broader and weaker than the Q branches which are char- 
acterised by narrow and strong absorption. The effects of Q 
branches in the spectra are thus often more apparent than 
the effects of P and R branches. This is seen in both figures 
[7| and [5] . The absorption features correspond to the places 
in which the Q branch lines are densest. Usually this corre- 
sponds to the position of the band centres. In some cases the 
line positions of a Q branch lines will increase in wavelength 
with J and at some point reverse and decrease in wavelength 
with J, in a similar fashion to a band head. For example the 
feature at about 3.58 fim in the synthetic spectrum is at- 
tributable to a Q branch of the (0,3,lf,l)— >(0,4,2e,l) band 
which displays such behaviour. 

This comparison of the Q branch line positions and in- 
tensities in figures |7| and |H| with the synthetic spectra al- 
lows us to attribute some of the absorption features in the 
observed spectra to Q branches. In figure the deep ab- 
sorption at 3.567 /xm in the spectrum of WZ Cas is at- 
tributable to the Q branches of the (0,0,0,0)^(0,1,1,1), 
(0,1,1,0)^(0,2,2,1), (0,2,2,0)->(0,3,3,l) bands and higher 
hotbands. The feature at 3.584 (jm is attributable in the 
most part to the Q branches of the (0,2,0,0)— ► (O, 3, 1,1) 
and the (0,3,1,0)^(0,4,2,0) bands. The line spacing of the 
(0,1,1,0)— ►(0,2,0,1) is large and consequently the intensity 
of the Q branch is spread over a larger wavelength range so 
is less detectable. Beyond 3.585 /im the increasing density 
of bands coupled with the errors on the ab initio data make 
it difficult to attribute ab sorption features to specific bands. 
The synthetic spectra of lAoki et al.l <ll998f) reproduce the 
3.56-3.63 /im band, but because of the limitations of their 
linelist they were unable to fully explain the various fea- 
tures in this region. The match of the synthetic spectra to 
the observed spectrum would be much improved if better 
frequency data was used. It is possible to substitute exper- 
imentally determined H CN/HNC energy lev els and hence 
line frequencies into the lHarris et al.l i2002tj) linelist. This 
would improve the line frequencies for all the lines corre- 
sponding to transitions between experimentally determined 
energy levels. However, this would require the manual as- 
signment of the "approximate" quantum numbers Wi, V2, «3 
and I to many of th e 125 000 energy levels includ ed in the 
lHarris et all ll2002bl) linelist. iBarber et all <2002h have al- 
ready assigned about a quarter of the levels we needed to 
be able to insert experimentally determined energy levels. 

The spectral region covered by figure |H| is slightly more 
complex, because there is significant absorption from both 
HCN and CS. The band heads of the CS 0->l, l-»2 and 
2^3 li e at 3.868 ym, 3.90 9 /mi and 3.9497 /im, respec- 
tively JChandra et al] |l99Sft . Three absorption features in 
the observed and synthetic spectrum of WZ Cas caused by 
these band heads are clearly visible and have previously 



been commented upon bv lAoki et al.l Jl998ft in WZ Cas 
and TX Psc and also bv lJorgensen et al.l J200(J) in TX Psc. 
As far as HCN absorption is concerned the 3.85 /mi fea- 
ture in the spectrum of WZ Cas is attributable to the Q 
branch of the (0,1,1,0)— >(1, 0,0,0) band. The absorption at 
about 3.88 microns can be confidently associated with the Q 
branches of the (0,2,2,0)^(1,1,1,0) and (0,4,2,0)^(1,3,3,0) 
bands. There is also a strong case for the 3.9 fim absorp- 
tion being caused by the unusually high density of Q branch 
lines in the (0,4,0,0)— ►(1,3,1,0) band. It is also likely that the 
(0,2,0,0)— ►(1,1,1,0) band is responsible for the weak absorp- 
tion at 3.857 fim in the spectrum of WZ Cas. The increasing 
density and decreasing strength of both HCN and CS lines 
coupled with the errors of the ab initio HCN data hinder 
further identifications longwards of 3.905 fim. 

In this work we have assumed that the spectrum of 
H 13 CN is identical to that of H 12 CN. In reality this is not 
the case, the vibrational ba nd centres and rot ational con- 
stants are slightly different. Maki et alJ ll200Cf) have mea- 
sured band centres and rotational constants for H 13 CN in 
the laboratory. These measurements found the band centres 
of the H 13 CN fundamentals to be lower in wavenumber than 
those of H 12 CN. The band centre of the fundamental of the 
C-N stretch mode differers by ~ 1.6% between the two iso- 
topomers, the bending mode by ~0.8% and the H-C stretch 
mode by ~0.5%. If the abundance ratio of 13 C to 12 C is 
high in a star, then one would find that the absoption bands 
of HCN would have slightly different shapes. Furthermore 
there would be additional Q branches of H 13 CN in the spec- 
tra, which would be slightly displace from those of H 12 CN. 
When corrected for the change in center of mass the H 13 CN 
electric dipole moment is identical to that of H 12 CN. Con- 
siquently the transition dipoles of the H 13 CN should be close 
to those of H 12 CN. To gain a good approxomation to the ab- 
sorption of H 13 CN the experimental energy levels of H 13 CN 
cou ld be subsitu t ed for the theoretical levels of H 12 CN in 
the lHarris et al.l j2002bl) linelist. This , however, must also 
wait until sufficent energy levels of the lHarris et alJ (|2002b) 
linelist are assigned. 

5.2.2 TX Psc. 

Figure El shows th e ISO SWS spectrum of TX Psc, taken by 
lAoki et alJ il998h . with our best fit synthetic spectra. One of 
the synthetic spectra is calculated with the full HNC/HCN 
opacity and the other is calculated with only HCN opacity. 
Both of the synthetic spectra are calculated with T B g = 3050 
K, log(C/0) = 0.03, V t = 3.5 km s" 1 and log(Z/Z Q )=-0.3. 

The T e ff and C/O ratio est i matio ns for TX Psc of 
lAoki et alJ (Il998h . lj0rgensen et all fcOOOl) and those of this 
work are summarised in table |3] together with interfero- 
metric estimates of T e fi ■ All the values of T e ff deduced 
from the fitting of synthetic spectra agree to within 100 
K, but the values obtained by interferome tric means are 
upto 200 K cooler. Jorgens en et alJ {2QQQ0 reported that 
their estimate of the T e fl of TX Psc changed from 3100 K 
to 3000 K between the two separate IS O observati o ns on e 
taken by themselves and one taken by lAoki et alJ <199Sft . 
However such temperature variations are difficult to con- 
firm as the error in T e g of our fit is around 100 K. The 
value of log(Nc/No) =0.015 from our work is closer to 
the lower log(A r c/A r o) =0.010 value of the two good fit 
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Figure 7. The observed spectrum of WZ Cas over the region of 3.55-3.61 (im, with a synthetic spectrum at T=2800 K, log(Nc / No) of 
0.003 and Vj* = 3.5 km s _1 . A constant is subtracted from the flux of the synthetic spect rum to stagger it from the observed spectrum. 
The Q branch lines in this spectral region in the ab initio data of lHarris et alj i2002b!) are plotted at the bottom of the figure with 
line length in proportion to intensity. The experimental band centres of the (0,0,0,0)— >(0, 1,1,1) band and its hot bands are indicated 
by the upper series of marks and digits. Likewise the lower series of marks and digits identify the band centres of the ab initio data of 
lHarris et alj l2002bl) . The data for these bands, together with its index are tabulated in table |21 



Table 2. Laboratory and ab initio band centres for bands with Ai>2 = l,Ai)3 = 1. 



index a 


(v[,v' 2 ,l',v' 3 ) 




E' c cm" 1 


E" cm- 1 


Theory ^ A c fim 


Exp c A c lira 


1 


(0,1,1,1) 


(0,0,0,0) 


2814.00 


0.0 


3.5537 


3.5643 


2 


(0,2,0,1) 


(0,1,1,0) 


3510.99 


715.89 


3.5777 


3.5841 


3 


(0,2,2,1) 


(0,1,1,0) 


3527.86 


715.89 


3.5562 


3.5652 


4 


(0,3,1,1) 


(0,2,0,0) 


4207.62 


1414.92 


3.5808 


3.5845 


5 


(0,3,1,1) 


(0,2,2,0) 


4207.62 


1432.20 


3.6031 


3.6040 


6 


(0,3,3,1) 


(0,2,2,0) 


4242.07 


1432.20 


3.5589 


3.5659 


7 


(0,1,1,2) 


(0,0,0,1) 


4891.95 


2100.58 


3.5825 


3.5953 


8 


(0,4,0,1) 


(0,3,1,0) 


4891.76 


2114.34 


3.6005 


3.6041 


9 


(0,4,2,1) 


(0,3,1,0) 


4907.37 


2114.34 


3.5803 


3.5849 


10 


(0,4,2,1) 


(0,3,3,0) 


4907.37 


2149.03 


3.6254 


3.6243 


11 


(0,4,4,1) 


(0,3,3,0) 


4957.36 


2149.03 


3.6123 




12 


(0,2,0,2) 


(0,1,1,1) 


5586.50 


2814.00 


3.6069 


3.6151 


13 


(0,2,2,2) 


(0,1,1,1) 


5602.92 


2814.00 


3.5856 


3.5965 



a This is the index of the band used in figure 171 



^ Calculated from the ab initio data of [ larris et al 
From the laboratory measurements of 



of IJeireensen et alj 


2000), than the higher values of 0.041- 


0.068 


lAoki et alllf 


!28j; I.Torgensen et, all2000 ; 





2000). In our synthetic spectra the HCN absorption features 
are very weak for values of log(Nc/No) greater than 0.02, 
see figure^] So we believe the lower log(Nc/No) values for 
TX Psc are closer to the real values. 

There are dramatic differences between the two syn- 
thetic spectra shown in figure |2l the broad 3.1 /im absorp- 
tion feature in the synthetic spectrum with HCN and HNC 



opacity is much weaker in the synthetic spectrum calculated 
with only HCN opacity. At T cff = 3000 K the model at- 
mospheres are very sensitive, slight changes in temperature 
and log(A r c/A r o) can have a dramatic effect on the syn- 
thetic spectra, for example see figure 0] and As discussed 
in section 15.11 a model without HNC opacity has a higher 
temperature and pressure at a given optical depth than a 
model with the full HCN/HNC opacity. Hence a reduction 
in opacity caused by neglecting HNC opacity can affect the 
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Figure 8. The observed spectrum of WZ Cas over the region of 3.8-4.0 /jm, with a synthetic spectrum at T=2800 K, \og{Nc / No) of 
0.003 and Vt = 3.5 km s" 1 . A constant is subtracted from the fl ux of the synthetic sp ectrum to stagger it from the observed spectrum. 
The Q branch lines in this spectral region in the ab initio data of lHarris et alJ i2002bl) are plotted at the bottom of the figure with line 
length in proportion to line intensity. The experimental band centres of the (0,1,1,0)— ►(1,0,0,0) band and its hot bands are indicated 
by the upper series of marks and digits. Likewise the lower series of marks and digits identify the band centres of the ab initio data of 
lHarris et alJ 12002 bl) . The downward pointing arrows indicate the position of CS band heads. The data for the HCN bands with their 
index is tabulated in table 151 
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Figure 9. Synthetic spectra of carbon rich atmospheres at T e fj =3050 K with a Vt=3.5 km s 1 and log(Nc/No) =0.015 in the region 
2.7 -4.0 lira. The spectra were calculated with the full HCN/HNC linelist of lHarris et"aT] 1 2002b ) and with just HCN lines selected from 
the lHarris et alJ 12002b!) linelist, see text. For comparison the observed spectrum of TX Psc is shown. We have also indicated with arrows 
the positions of CS and CH band heads. The resolution of these spectra is 1.6 nm. 
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Table 3. Laboratory and ab initio band centres for bands with A^i = l,Av2 = — 1. 



index a 




{v'{,v'^,l" ,v'^) 


E' c cm 1 


E'J cm" 1 


Theory ^ A c fim 


Exp c A c (im 


1 


(1,0,0,0) 


(0,1,1,0) 


3307.70 


0.0 


3.8588 


3.8469 


2 


(1,1,1,0) 


(0,2,0,0) 


4004.28 


1414.92 


3.8620 


3.8569 


3 


(1,1,1,0) 


(0,2,2,0) 


4004.28 


1432.20 


3.8879 


3.8795 


4 


(1,2,0,0) 


(0,3,1,0) 


4686.28 


2114.34 


3.8881 


3.8898 


5 


(1,2,2,0) 


(0,3,1,0) 


4702.19 


2114.34 


3.8642 


3.8673 


6 


(1,2,2,0) 


(0,3,3,0) 


4702.19 


2149.03 


3.9167 


3.9132 


7 


(1,3,1,0) 


(0,4,0,0) 


5366.17 


2801.46 


3.8991 


3.9003 


8 


(1,0,0,1) 


(0,1,1,1) 


5394.43 


2814.00 


3.8753 


3.8638 


9 


(1,3,1,0) 


(0,4,2,0) 


5366.17 


2817.27 


3.9233 


3.9236 


10 


(1,3,3,0) 


(0,4,2,0) 


5400.32 


2817.27 


3.8714 


3.8782 


11 


(1,3,3,0) 


(0,4,4,0) 


5400.32 


2867.21 


3.9477 


3.9480 



This is the index of the band used in figure 151 



Calculated from the ab initio data of Harris ct al. 



From the laboratory measurements of Ma ki et alj 11996 



Table 4. Effective temperatures, log(A r c /No) and surface gravities for TX Psc. 



Work 


T e ff 


log(N c /N ) 


log(g) 


Notes 


Lambert et al. ( 1986) 


3030 


0.012 




T e ff from colour index, log(7Vc/./Vo) from models. 


Ouirrenbach et al. (19941 


2805±126 






intcrfcromctry 


Dvck et al. ( 19961 


2921±60 






interferometry 


Aoki et al. (19981 


3100 


0.041 


0.0 


Recalibrated IRFM determination of Ohnaka & Tsuii 


Joreensen et al. (2000") 


3100 


0.01 


-0.5 


model 


Jore-ensen et al. (20001 


3100 


0.041 


0.0 


model 


j0rsrensen et al. (2000) 


3000 


0.01 


-0.5 


model 


Joreensen et al. (2000) 


3000 


0.041 


0.0 


model 


Ohnaka et al. (2000) 


3000 


0.03 


0.0 


model 


Ohnaka et al. (20001 


3100 


0.068 


0.0 


model 


Ohnaka et al. (20001 


3000 


0.041 


-0.5 


model 


This work 


3050 


0.015 


0.0 


model 



structure of the atmosphere to the extent at which HCN ab- 
sorption is dramatically reduced. As in the case of WZ Cas 
both HCN and HNC are required to explain the observed 
spectra. Marked on figure |5] are the 0— *T, 1— >2 and 2^3 
CS and CH band heads, which were identified previously by 
lAoki et alJ Jl998f) . 

Although the synthetic spectra calculated with 
HCN/HNC opacity fits the observed spectrum well between 
2.9 to 4.0 fim, there is a slight under-prediction of the 
strength of the absorption between 3.0 and 3.2 fim. Between 
2.7 and 2.9 fim the flux from the synthetic spectrum dips 
to below the spectrum of TX Psc. The absorption across 
this region is effected by the tail of the CO Av = 2 bands. 
The effect of this CO band is evident in the synthetic spec- 
tra shown in figure This would imply that there is too 
much CO within our model atmosphere. We have looked 
at other potential ways of reducing the abundance of CO, 
including changing the surface gravity and increasing and 
decreasing the value of log(iV)sr/iVo) above and below the so- 
lar value. It is clear that both these parameters do affect the 
HCN absorption features in the final synthetic spectrum, but 
we have not found a better fit. The elemental abundances 
within carbon stars are not well known, so \og(Nc/Nn), 
log(A r N /A r H), log(N /Nn) or indeed log(Ar H e/Ar H ) are not 
well constrained by observational data. With so many free 



parameters it is possible there are other good fits that we 
have not found. 



6 CONCLUSION. 

We have calculated new model atmospheres and syn- 
thetic spectra for car bon stars, us i ng the recent ab initio 
HCN/HNC linelist of lHarris et ail i2002tJ> . However selec- 
tively removing HNC lines from the HCN/HNC linelist and 
recalculating our models, we have shown that HNC opacity 
can have a small but significant effect upon the structure of 
carbon star atmospheres. The general shape of the pressure 
temperature curves of our model at mospheres w h ich ex clude 
HNC are consistent w ith those of lAoki et alJ lll998ll and 
Ij0rgensen et alj i2000ll . However our models with both HCN 
and HNC will differ f rom those of lAoki et ail <199Sh and 
Ijelrgensen et alj i2000ll . especially at low effective tempera- 
tures (T e ff< 2900 K) and low \og(N c /N ) (\og(N c /N ) < 
0.01). 

We have obtained good fits of our synthetic spectra to 
the observed spectra of WZ Cas and TX Psc. These fits 
have provided a new independent estimate of the C/O ra- 
tio and T e ff of these stars. In general our estimate of the 
T c ff and C/O ratio for WZ Cas and TX Psc are in agree- 
ment with earlier works. For WZ Cas our T e ft and C/O ratio 
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are both slightly lower than earlier estimates, this is consis- 
tent with an increase in overall opacity within the model due 
to the inclusion of HNC. For TX Psc our C/O ratio is to- 
ward the lower end of the earlier estimates of the C/O ratio. 
Our models using values of \og(Nc/No) greater than 0.02 
favoured by earlier works cannot reproduce the observed 
spectra of TX Psc. However by changing the many free 
parameters such as log g \og(Nc/No) , log(Nj$ /No) and 
T c ff i t may be possib l e to o btain other good fits, for ex- 
ample Ijprgensen et alj JgOOO) obtained two good fits with 
different values of \og(Nc/No) , T e fi and log g. 

By comparing our synthetic spectra calculated with and 
without HNC opacity, we have identified the contribution of 
HNC at 2.9 /xm to the spectrum of WZ Cas for the first 
time. Absorption resulting from the Q branch lines of the 
HCN Av 2 = l,Av 3 = 1 and Avi = l,Av 2 = -1 bands 
has been identified in the spectrum of WZ Cas. It has been 
shown that to correctly model the atmospheres of C-stars, 
opacity arising from HNC must be accounted for. 

It may be possible to co rrect for systematic e rrors in line 
frequencies of the linelist of lHarris et alJ (1200 2b) by substi- 
tuting experimentally determined energy levels for their ab 
initio counterparts. This would give very accurate line fre- 
quencies for the lines corresponding to transitions between 
the experimentally determined energy levels. Similarly by 
substituting experimentally determined H 13 CN energy lev- 
els for their ab initio H 12 CN counterparts, it should be possi- 
ble to better account for the spectrum of H 13 CN. This, how- 
ever, must await the assignment of approxim ate quantum 
number s to the energy levels of the linelist of lHarris et alJ 
J2002bl) . 
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